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Studying Yang-Mills field and gravitational field in class A Bianchi spacetimes, we find that 
chaotic behavior appears in the late phase (the asymptotic future). In this phase, the Yang-Mills 
field behaves as that in Minkowski spacetime, in which we can understand it by a potential picture, 
except for the types VIII and IX. At the same time, in the initial phase (near the initial singularity), 
we numerically find that the behavior seems to approach the Kasner solution. However, we show 
that the Kasner circle is unstable and the Kasner solution is not an attractor. From an analysis 
of stability and numerical simulation, we find a Mixmaster-like behavior in Bianchi I spacetime. 
Although this result may provide a counterexample to the BKL (Belinskii, Khalatnikov and Lifshitz) 
conjecture, we show that the BKL conjecture is still valid in Bianchi IX spacetime. We also analyze 
a multiplicative effect of two types of chaos, that is, chaos with the Yang-Mills field and that in 
vacuum Bianchi IX spacetime. Two types of chaos seem to coexist in the initial phase. However, 
the effect due to the Yang-Mills field is much smaller than that of the curvature term. 

I. INTRODUCTION 

■ The Yang-Mills field is one of the central objects in particle physics. When we discuss a unification of fundamental 
interactions, non-Abelian gauge fields, which are a generalization of the Yang-Mills field, play a key role. Hence we 
should study more carefully the Yang-Mills field in cosmology, in particular in the early stage of the universe. If 
the field is random, it is described by "radiation", which is simple to deal with. The Yang-Mills field may also be 
homogeneously distributed by some symmetry-breaking mechanism. So far, there has not been much study of the 
' homogeneous Yang-Mills field in the context of cosmology. In fact, the homogeneous source- free Yang- Mills field has 
been studied in the hope that a non-perturbative treatment may allow a better understanding of the vacuum state 
of the Yang-Mills field, despite the fact that strong and weak interactions have no classical counterpart. It is also 
r~ p interesting to study the Yang-Mills field from the viewpoint of a dynamical system because it shows chaotic behaviors 
in Minkowski spacetime 0, 0, IE El HI E| ( see also a good review 0). The homogeneous Yang-Mills field evolves as a 
^Jy particle in some potential. A particle in such a potential (see Fig. |2Jfor the potential shape) is known to show chaotic 
behavior by the study of a billiard system. The Yang-Mills field, however, is integrable in the Friedmann universe 0. 
This result suggests that symmetry of the universe spoils its chaotic property because a number of dynamical degrees 
of freedom decreases due to the assumed symmetry [ij. In axisymmetric Bianchi I spacetime, the Yang-Mills field 
again shows chaotic behavior lOHTj. Which spacetime shows the chaotic behavior of the Yang-Mills field? Is chaos 
a universal phenomenon in generic spacetimes? With these questions in mind, we search the Einstein- Yang-Mills 
(EYM) system for more generic Bianchi type spacetimes [T^ . 

We should also stress another important point in the present system. We know that the past asymptotic behavior 
toward the initial singularity of vacuum Bianchi VIII and IX spacetimes is chaotic [Fj LLi UM UH LLlI ■ They are called 
the Mixmaster universe [f3 |. The other type Bianchi spacetimes with a perfect fluid do not show chaotic behavior. For 
example, a Bianchi I spacetime with a perfect fluid approaches monotonously the Kasner spacetime, which is a shear 
dominant universe. The inclusion of a magnetic field, however, makes the dynamics very complicated. Bianchi I, II, 
VIn spacetimes with a perfect fluid and a homogeneous source-free magnetic field behave like a Mixmaster universe 
[l8l fl9L |20| . We should then study whether the past asymptotic behavior of Bianchi spacetimes with the Yang-Mills 
field shows similar chaotic behavior. 

We also mention the BKL (Belinskii, Khalatnikov and Lifshitz) conjecture, which predicts a fate of singularities 
in general inhomogeneous spacetimes |2l|,|23|. The BKL conjecture naively says that any spatially different points 
decouple each other near initial singularity. Strictly speaking, it consists of the following three statements: 
First, a spacetime near singularity becomes spatially homogeneous at each spacetime point, which can be described 
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explicitly by using a suitable local coordinate system. 

Secondly, the contribution of a matter to the system becomes negligible near singularity compared with those of 
shear and curvature terms, although the matter energy density may blow up into initial singularity. This asymptotic 
spacetime is given by a vacuum Bianchi model. This is why vacuum models are often studied. 

Thirdly, the most generic spatially homogeneous models near singularity is Bianchi IX spacetime. These statements, 
then, yields that a spacetime near initial singularity should be describe d by the Mixmaster universe. Now let us 
look back at the case of Bianchi I, II, VIo model with magnetic field |18l fl9L l2Cj. In their cases, the energy density 
of magnetic field, which is normalized by a Hubble expansion parameter, does not vanish toward a singularity. It 
sometimes shows a spiky growth like a curvature term in the Mixmaster universe. This breaks the second statement 
of the BKL conjecture. Hence, such a model may be regarded as a counterexample to the BKL conjecture. It is 
interesting to study whether Bianchi I (or more generically Bianchi IX) spacetime with the Yang-Mills field also give 
a counterexample to the BKL conjecture. 

A natural question arises in the case of Bianchi IX spacetime: Does a combination of two types of chaos, that is 
chaos with the Yang-Mills field and chaos in the Bianchi IX spacetime, strengthen each chaotic behavior? We now 
proceed to analyze the Yang-Mills field in the Bianchi IX spacetime. 

Finally, it should be noted that it is difficult to discuss chaos in general relativity. Many standard indicators 
of chaos, such as the Lyapunov exponent, depend on the choice of coordinate systems. In order to define chaos 
independently from the coordinate choice, many efforts have been done so far. A fractal basin approach is well known 
for the coordinate independent indicator of chaos llo . Il7l . Using the Lyapunov exponent, recently, a coordinate 
independent definition of chaos has also been proposed |23|. In this paper, we adopt a fractal basin approach to 
determine whether a system is chaotic or not. 

We use the units of c = h = 1. 



II. BASIC EQUATIONS 



The action of the EYM system is 

S = 



-9 



R 



1 p(A) p(A)nv 



(2.1) 



where R is the scalar curvature of the metric g^ v and F^v ' is the field strength of Yang-Mills field. " A v describes 
the components of the internal space. The variation of the action with respect to a metric g^ v and a vector potential 



gives two basic equations, which are the Einstein equations 

Gfj, u = 87rGTp„, 
1 



■9ym 



and the Yang-Mills equations 



F"".„ - A v x F^ = , 



(2.2) 
(2.3) 

(2.4) 
a\^ta with 



where <7ym is the self-coupling constant of the Yang-Mills field, and where F^" — F^'^ta and A u 
ta being the SU(2) basis. When we set our units as 87tG/<7 Y m = 1 as we U as c = Ti = 1, the basic equations become 
free from the value of gvM- 

We adopt the orthonormal- frame formalism, which has been developed in j24[. In Bianchi- type spacetimes, there 
exists a 3-dimensional homogeneous spacelike hypersurface S t that is parameterized by a time coordinate t. The 
timclike basis is given by eo = df . The triad basis {e } on the hypersurface S t is defined by the commutation 
function 7 c Qb as 



[e a , e b ] = i c ab e c . 



(2.5) 



It is convenient to decompose 7 c ab into the geometric variables (the Hubble expansion H, the acceleration it a , the 
shear a ab , and the vorticity uj ab ), and the variables tt a (the rotation of e a ), and the variables a a and n ab , which 
distinguish the type of Bianchi models, i.e., 



7 ob 
7° 0a 



H5 b a - eV(c c + n c ) 



-2e„ 



l c ab = tabdn dc + a a S b c - a b S a c 



(2.6) 
(2.7) 
(2.8) 
(2.9) 
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In Bianchi spacetimes, the vorticity u a b always vanishes because they are hypersurface-orthogonal. 

In order to analyze the EYM system, we have to write down the basic equations. For the generic situation, however, 
we have so many variables. Hence we shall study the simplest case. First, using the gauge freedom of the Yang-Mills 



field, we set Aq \t) = 0, which simplifies the vector potential as A — A)^' (t)TA<jJ a , where ui a is the dual basis of e a 



(A), 



In the class A Bianchi type, a a vanishes. We can also diagonalize n a b, i.e. n a b — diagjni, 712, 713}, using the 
remaining freedoms of a time-dependent rotation of the triad basis. Then, we can show that if a a b, A^ , and A^ do 
not have initially off-diagonal components, the equations of motion guarantee that those variables always stay diagonal 
for class A Bianchi spacetime (see Appendix). As a result, tt a vanishes, and the number of basic variables reduces 

from 21 (12 for spacetime [H, N a , a a b, fi a ] and 9 for the Yang-Mills field [^i^]) to 10 (7 for spacetime [H, N a , a aa ] 
and 3 for the Yang-Mills field [Y^]). This is not valid for the class B Bianchi spacetimes. Therefore, in what follows, 
we discuss only the class A Bianchi spacetimes. 

In order to discuss the dynamics, it is also convenient to introduce the Hubble normalized variables, which are 
defined as follows: 



^ab — 



TV = — 

H 



We also introduce the conventional shear variables of the Bianchi spacetime as 



= 5 ( S 22 



£33), £- 
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The diagonal components of the Yang-Mills field potential are described by new variables a,b,c as 

4* 
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(2.10) 



(2.11) 



(2.12) 



A new time variable r defined by dr — Hdt is introduced, r denotes the e-folding number of the scale length. We 
choose the origin of r when H = 1. 

Here we state the basic equations explicitly. They consist of the generalized Friedmann equation, the dynamical 
Einstein equations, and the Yang-Mills equations. The generalized Friedmann equation, which is the Hamiltonian 
constraint equation, is described as 



CVm + K = 1 



(2.13) 



where JIym is the density parameter of the Yang-Mills field, i.e. the Hubble normalized energy density of the Yang- 
Mills field, which is given by 



YM 
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and where K is the Hubble normalized curvature, which is defined by 

( 3 )i? 1 



K 



{Nf 



2{N 1 N 2 + N 2 N z + N i N 1 )} 



(2.14) 



(2.15) 
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Note that the positive 3-curvature corresponds to K < 0. Then, except for the Bianchi type IX, K > 0. The energy 
density of the Yang-Mills field is always positive definite. Thus we find that E 2 , fiyivb and K are restricted as 
< E 2 < 1, < J1 YM < 1, and < K < 1, except for type IX. 
The dynamical Einstein equations are 

(2.16) 
(2.17) 
(2.18) 
(2.19) 
(2.20) 
(2.21) 
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where 



S+ = ~{{N 2 ~N 3 ) 2 - Ni(2Ni - N 2 - N 3 )} 



S- 



2\/3 



(N 3 - N 2 )(N 1 - N 2 - N 3 ) , 



which depend on only N a , and 



n 4 



n 



2 {a' + (-253+ + l)af + {&' + (E+ + V3E_ + 1)&} + [c' + (E+ - V3E_ + l)c}" 

2 („ 10 + |) 2 + (w2t+i) V( W3C+ g) ! 

{&' + (E+ + V3E_ + l)fe} 2 + |c' + (53+ - V3E_ + l)cV 
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which are the Hubble normalized anisotropic pressures of the Yang-Mills field. 
The Yang-Mills field equations are 

a" = (Y?-K-l)a' 

+ (E 2 -5K + AY? + - 4E+ - Ni(N 2 + N 3 ) + 2U + )a 

_1(7V 1 + N2 + N3 ) bc - _L a (62 + c 2 } 

b" = (E 2 - K-\)b' 

+ (E 2 ~5K+ (E+ + \/3E_)(E + + \/3E_ + 2) - 7V 2 (iV3 + Ni) -11+ - y/S£L.)b 

- — {N x +N 2 + N 3 )ca - j^b(c 

c" = (E 2 -X-l)c' 

+ (E 2 -5K+ (E+ - V3E_)(E+ - V3E_ + 2) - N 3 (Ni + N 2 ) -U+ + V3U^)c 

~(N! +N 2 + N 3 )ab - j^c{a 2 + b 2 ). 



2 ' a 2 ) 



(2.22) 
(2.23) 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



The basic variables are {H, E + , E_, Ni, N 2 , N 3 , a, b, c}, whose equations are given by Eqs. I|2.16|l . I|2.17|l . 
QHfy . QJB , ijOljl . f2~2"TJ) . (|2~2^jl . l(2~27| . J22HI, with one constraint equation fmty . 

For all class A Bianchi spacetimes except for type IX, we numerically find two typical phases in the evolution of 
the universe; i.e., an initial phase near a big-bang singularity, and a late phase, when a potential picture for the 
Yang-Mills field becomes valid. We will show typical behaviors in each phase in Fig. including a transient phase. 
In what follows, we discuss the details of each phase in order. 



III. THE LATE PHASE 



In an axisymmetric Bianchi I model, Darian and Kunzle show that chaos appears in the dynamics of the Yang-Mills 
field [10|. They adopt a potential picture to understand its chaotic behaviors. Can we adopt the some approach to 
analyze the present system? To compare our analysis with their results, we shall transform our tetrad basis to the 
Cartesian coordinate basis. The Yang-Mills field potential a, b, c should be replaced with A, B, C as 

A = £a, B = £b, C = £c, (3.1) 

where I = £(t) is a scale factor defined by H = £/£ and given explicitly as £(t) — £$e T ~ e T . £q is the initial value of 
the scale factor at r = 0. Without a loss of generality, we can set £q — 1 by rescaling because the Hubble-normalized 
curvature K does not depend on £{t). 

When we discuss the dynamics of the Yang-Mills field, adopting a potential picture makes it convenient to see 
the existence of chaos, just as in Darian and Kunzle. For this purpose, we change our time coordinate from r to 
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FIG. 1: Two typical phases in the Bianchi I spacetime: (a) one is the initial phase (E-dominance) and (b) the other is the 
late phase (f^YM-dominance) , where E, Hym, and K denote the Hubble normalized shear strength, the Hubble normalized 
energy density of Yang-Mills field, and the Hubble normalized curvature, respectively. These are constrained by the generalized 
Friedmann equation: !!ym + E 2 + K = 1 . 

a conformal time 77, which is defined by dr\ — £(t)~ 1 dt = H~ 1 e~ T dr. Using this new coordinate, the Yang-Mills 
equations are written as 



d 2 A 
drj 2 

d 2 B 
drj 1 

d 2 C 
drj 2 



H 2 e 2T (-4K + AY? + - 4S+ - Ni(N 2 + N 3 ) + 2U+)A 

dV 

-He T {N 1 +N 2 +N 3 )BC-—, (3.2) 

H 2 e 2r (-4K + (S+ + V3S_)(E+ + V3£_ + 2) - N 2 (N 3 + Ni) -11+ - V3U-)B 

dV 

-He T (Ni + N 2 + N 3 )CA - — (3.3) 

H 2 e 2T (-4K + (S+ - V3£_)(£ + - V3S_ + 2) - N 3 (Ni + N 2 ) + V3n_)C 

dV 

-He T (N 1 +N 2 + N 3 )AB- — , (3.4) 



where a potential V is defined by 



dC 



V(A, B, C) = ^(A 2 B 2 + B 2 C 2 + C 2 A 2 ) , (3.5) 



whose shape is schematicallyshown in Fig. [21 It is known that a test particle moving in such a potential shows a 
chaotic behavior 0, 0, 0> Then, if we can adopt the potential picture, that is, if the potential term becomes 

dominant, we expect chaos in the dynamics of the Yang-Mills field. 

To show it, we divide the r.h.s. of Eqs. I|3.2|l . I|3.3|) . and (|3.4ll into three terms. The first one is proportional to the 
field itself (A, B, or C) and coupled to the background geometrical variables. The second one gives a coupling with 
other components of the Yang-Mills field through N a (the variables classifying Bianchi types). The third one contains 
only the Yang-Mills field variables. We then call these three terms an oscillation term, a Bianchi-mixing term, and a 
potential term, respectively. 

If we can ignore the first and second terms, we can discuss the dynamics only via the potential V and adopt a 
potential picture for the present dynamics. In what follows, analyzing the above three terms, for each Bianchi model, 
we show that a potential term becomes dominant and a potential picture is valid in the late phase for some Bianchi 
models. 



A. Bianchi type I, II, and VIo spacetimes 



We first show that a potential picture is valid in the Bianchi I, II, and VIo spacetimes in the late phase by evaluating 
each term in Eqs. $£2$ , JjHjjl, and J32J. 
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• Since K > 0, E 2 > 0, and flyu > 0, from Eq. (|2.13() . we find that JIym, E 2 , and K are bounded as < Hym < 1, 
< E 2 < 1, and < K < 1. E ± is also restricted as -1 < E± < 1. 

• Eq. (|2.16(l guarantees that H damps exponentially. Its exponent index is in the range of —1 (curvature 
dominance) to —3 (shear dominance). 

• II± is bounded as follows: 
Using JJym, n + is rewritten as 



n + — -Qym + 2 



{a' + (-2E+ + l)a} 2 



Nia- 



(3.6) 



which gives the constraint of — 1 < ±1 + . n + is also rewritten as 



IE 



2Q 



YM — 7T 



{&' + (E+ + VS£. 
+{c' + (E+ - V3E_ - 



+ 1)&} 2 
l)c} 2 4 



^26 + -^ca 



N 3 c- 



H 



-ab 



(3.7) 



which gives another constraint of II + as II + < 2. As a result, II + is bounded as — 1 < II + < 2. II_ is similarly 
bounded as -a/3 < n_ < V3. 

• In the case of Bianchi types I, II, and VIo , we can show that N a is bounded, that is, 

— In Bianchi I spacetime, N a = 0. 

— In Bianchi II spacetime, K = yj-ZV 2 . Hence iV a is bounded as — 2\/3 < Ni < 2a/3. 

— In Bianchi VIo spacetime, K — yj (^Vi — N2) 2 , and Ni and N% have different signs, i.e., NiN% < 0. Hence, 
N a is bounded as -2\/3 < Ni,N 2 < 2x/3. 

From the above discussions, we find that both the oscillation and the Bianchi-mixing terms in Eq. (|3.2|l damp and 
then that the potential term dominates with time, if the Yang-Mills field A, B and C do not decay exponentially. 
There is one exceptional case, in which H decays exactly as e~ T . Hence, except for such a special case, the Yang-Mills 
field behaves asymptotically (r — > 00) as a particle moving in the potential V(A, B, C). This potential shape is the 
same as that in Minkowski spacetime. Then, the Yang-Mills field will show a chaotic behavior just as in Minkowski 
spacetime 0. 

In order to confirm the above results, we also analyze the present system numerically. To show its chaotic property, 
we show its sensitive dependence to initial data. The potential (|3.5|l has three channels along the A, B, and C axes. 
We then set a window at a sufficient distance in each channel. When the orbit passes through one of these three 
windows, we assign three colors to the initial data in the following way: "bright gray," if the orbit passes through 
a window along A axis; "dark gray," if it falls down in a window along the B axis; and "black" for a window along 
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(a) (b) 

FIG. 3: Fractal basin boundaries of initial data in Bianchi I spacetime. We set three windows at \A\ = 2.8 (bright gray), |B| = 2.8 
(dark gray), and |C| = 2.8 (black) in three channels along A, B, and C axes, where the fates of final state of orbits are determined. 
Initial conditions are as follows: r = 0, H = 1, £+ = 0.02, E_ = 0.001, JVi = N 2 = N 3 = 0, a' = 0.3, b' = 0.5, d = 0.7. 
Further, 

(a) A(= a) and C(= c) vary from 0.6 to 1.0 . 

(b) An enlargement map of a part of Fig. El (a) [0.73 < A(= a) < 0.75 and 0.75 < C(= c) < 0.77]. 
B(= b) is obtained from Eq. (12.1311 . In each cases, a grid is 200 x 200. 



the C axis. If numerical simulation is broken by increasing numerical errors before the orbit reaches any window, we 
assign "white" to such initial data. We show our results in the Fig. |3] Fig. |3] shows that basin boundaries are fractal 
and orbits are sensitive to initial data, that is, orbits are chaotic. We find a relation between fractal basin boundaries 
and chaos similar to the other dynamical systems 0, 0] . 

In order to see a strength of chaos, we calculate a fractal dimension of the basin boundaries of Fig. |3 We adopt a 
box-counting method to evaluate the dimension. We find that the fractal dimension of the "black" region in Fig. |3| 
(a) is D = 1.687 (see Fig. |J. 

4.5 
4 
3.5 

log 10 « L 3 
2.5 
2 
1.5 

-0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 
log 10 i 

FIG. 4: The fractal dimension of the "black" region in Fig. El (a). The horizontal axis shows the logarithmic size of boxes 
(log 10 L), where L is normalized by a grid size, while the vertical axis gives the logarithmic number of boxes (log 10 til,) which 
contain black points. This figure shows the fractal dimension is D — 1.687. 
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B. Bianchi Type VIIo and VIII spacetimes 

In Bianchi VIIo and VIII spacetimes, it may be hard to show whether a potential picture is valid or not, because 
N a is not bounded. We may, however, be able to verify it by evaluating the "energy" E numerically: 

E= \ A ^ + l B 'i + + \( A2B2 + B2 ° 2 + C2A ^- (3 - 8) 

If the potential picture is valid, E will be asymptotically constant. We could successfully adopt the potential picture. 
We shall then evaluate E in the case of Bianchi VIIo and VIII spacetimes. 

Fig. (a) shows that E approaches a constant for a Bianchi VIo model, which was shown to be chaotic in the 
previous subsection. In Bianchi VIIo spacetime, E seems to approach a constant; therefore a potential picture may 
be valid, even if N a diverges (see Fig. 0(b)). In Bianchi VIII spacetime, however, if two of N a diverge, then E does 
not converge, and a potential picture may be no longer valid. However, we also find that when 7V does not diverge, 
E seems to approach constant, as in the case of Bianchi VIIo spacetime. It seems to depend on initial conditions 
whether N a diverges or not. We could conclude that in Bianchi VIII, chaos via a Yang-Mills potential appears for 
some initial conditions (i.e., iV a is bounded), but a different type of dynamical behavior is found for other initial 
conditions (i.e. N a diverges). We are not sure whether the latter case is chaotic or not. 
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(a) (b) (c) 

FIG. 5: Time evolution of E: (a) Bianchi VIo spacetime. E approaches constant for this chaotic system, in which a potential 
picture is valid, (b) Bianchi VIIo spacetime. E seems to converge within 20% fluctuations, (c) Bianchi VIII spacetime. E 
seems not to converge. 

We also analyze the sensitive dependence on initial data for Bianchi VIIo and VIII spacetimes. From Fig. we find 
that basin boundaries are fractal and the box-counting dimension of the " black" region in Fig. (a) is D = 1.714. 
We conclude that the EYM system in the Bianchi VIIo system is chaotic as expected from the evaluation of E. 

On the other hand, Fig. (for Bianchi VIII) looks different from the previous fractal structures. Most parts of 
Fig. 0are black. The box-counting dimension of the "black" region in Fig. 0(a) is D = 1.921. The box-counting 
dimension is much closer to D = 2 compared with the cases of Bianchi I and VIIo . Then Fig. seems not to be 
fractal. It seems that the oscillation and the Bianchi-mixing terms enlarge C, and then most orbits get into the C 
window. It may depend on the values of N a which window is chosen by most orbits. Figs. and support our 
expectation which is obtained by asymptotic behavior of E given in Fig. 

C. Bianchi type IX spacetime 

In the case of Bianchi IX, the behavior is very similar to Bianchi VIII. The most important difference is the 
asymptotic behavior of the spacetime. Since Bianchi IX spacetime is closed, once the Hubble normalized variables 
diverge, the Hubble parameter H eventually vanishes within a finite time, and the universe evolves into a big crunch. 
In that case, we cannot see anything about chaos in this period. 

Let us take a more detailed look. If curvature term K is negative, H vanishes rapidly, so that Yang-Mills field 
A, B, C do not have enough time to oscillate many times. We cannot determine whether such a system is chaotic 
or not. If one of N a is much smaller than the others, e.g., N3 <§; N\,N2, however, K becomes positive as K w 
(Ni — N 2 ) 2 /12, and then H may not vanish so rapidly, so that the Yang-Mills field A, B, C could have enough time 
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(a) (b) 

FIG. 6: Fractal basin boundaries of initial data in Bianchi VIIo spacetime. We set three windows at \A\ = 2.8 (bright gray), 
\B\ = 2.8 (dark gray), and |C| = 2.8 (black) in three channels along A, B, and C axes, where the fates of final state of orbits 
are determined. Initial conditions are as follows: r = 0, H = 1, £+ = 0.02, E_ = 0.001, Ni = 0.2, N 2 = 0.3, iV 3 =0, a' = 
0.3, b' = 0.5, d = 0.7 . Further, 

(a) A(= a) and C(= c) vary from 0.6 to 1.0 . 

(b) An enlargement map of a part of Fig. (a) [0.88 < A{= a) < 0.9 and 0.93 < C(= c) < 0.95]. 
B(= b) is given from Eq. 12.131 . In each cases, a grid is 200 x 200. 
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FIG. 7: Fractal basin boundaries of initial data in Bianchi VIII spacetime. We set three windows at \A\ = 2.8 (bright gray), 
\B\ — 2.8 (dark gray), and |C| = 2.8 (black) in three channels along A, B, and C axes, where the fates of final state of orbits 
are determined. Initial conditions are the following: r = 0, H = 1, £+ = 0.02, E_ = 0.001, Ai = 0.2, N 2 = 0.3, 7V3 = 
-0.2, a' = 0.3, b' = 0.5, c' = 0.7. Further, 

(a) A(= a) and C(= c) vary from 0.6 to 1.0 . 

(b) An enlargement map of a part of Fig. (a) [0.88 < A{= a) < 0.9 and 0.93 < C(= c) < 0.95]. 
B(= b) is given from Eq. 12.131 . In each cases, a grid is 200 x 200. 

to oscillate many times. If one of N a is dominant, N a is bounded as N a < 2\/3, using the generalized Friedmann 
equation l|2.13|l as long as K is positive. In this case, E seems to approach a constant value, just as in the case of 
the Bianchi VIII spacetime, so that a potential picture holds. Although the spacetime eventually recollapses into a 
big crunch, we can find a temporal chaos for some period. However, depending on its initial conditions, we may find 
divergence of N a , and then a negative K in the evolution of the spacetime. In this case, we will find a big crunch 
before the appearance of chaos in the Yang-Mills field. 

When a potential picture is valid, the oscillation term and the Bianchi term turn out to be negligible. In this case, 
there is not much difference between Bianchi I and IX spacetimes. We conclude that there is no enhancement effect 
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by spacetime chaos even in Bianchi IX spacetime because the vector potential behaves as in Minkowski spacetime. 



IV. THE INITIAL PHASE 



A. Bianchi type I spacetime 



From Fig. ^ wc nn d that the energy density of the Yang-Mills field (Uym) is damping to zero in the past direction 
near the initial singularity. The Bianchi I EYM system seems to approach the vacuum Bianchi I system, i.e., the 
Kasner solution. The Kasner solution is described by E 2 = 1 setting Uym = K = 0, which is called the Kasner circle. 
The Kasner circle is a past attractor in the case of Bianchi I spacetime with a non-tilted perfect fluid < 7 < 2 [24| , 
where 7 is a parameter of an equation of state p — (7 — l)/i. However, we have to be careful in drawing conclusions 
from our numerical analysis. In the case of the Bianchi IX model, the evolution is described by a sequence of Bianchi 
I, i.e. the spacetime first approaches to some Bianchi I spacetime but it bounces via a Bianchi IX potential back 
to another Bianchi I spacetime. This process repeats infinite times to the initial singularity. Even in the Bianchi I 
model, if we include a magnetic field, it is shown that the Kasner circle has a two-dimensional unstable manifold in 
the past direction, just as in the Bianchi IX model |19|. So it is important to study whether the Kasner circle is a 
past attractor or not in our Bianchi I EYM system. 

First, we consider the case of the Bianchi I model. From Fig. ^ one ma Y anticipate that the Kasner circle is a past 
attractor. Suppose that the spacetime approaches to one Kasner solution, that is K = 0, S^ym ~ 0, ^ — > ip (some 
constant), and E — * 1 (the Kasner circle), where E + = Ecos'I', E_ = Esin\I/. Near the Kasner circle, we obtain the 
following equation 



CI 



1 



YM 



{a' - (2Ecos* - l)a} 2 + U' - ^2Eco& 



cos I ^ + -TT I - I I b 



Nia + ^bc\ + (N 2 b + ^ca\ + (n 3 c + ^ab 



2Ecos f * — -7T 1 I ] r 



(4.1) 



which yields 



a' = (2Scos* - l)a, b' = 



2 Y.Cs j * + -TT ) - 1 > b, 



2EcOS I * - -TT 1 - 1 } r. 



(4.2) 



These equations lead to a situation in which at least one variable of a,b,c diverges exponentially in the limit of Bianchi 
I spacetime (E = 1 and * = ip) and one variable at least decays exponentially in the past direction. Using Eq. (|4.2|l . 
we evaluate the last terms in Eq. (|4.1|) . 



be 
H 
ca 
~H 
ab 
~H 



— oc exp [{2 cos(f/; + tt) + l}r] 
cx exp 
oc exp 



2 cos ( ip — -tt J + 1 f T 



2 cos [ ip + -tt ) + 1 } T 



(4.3) 



If 



< v < 



if t^tt < ip < |7r, and if — ^tt < ip < tt, then ab/H, bc/H, and ca/H damps to the past (t 



-00), 



respectively. As a result, ab/H, bc/H, and ca/H do not damp simultaneously, and then J]ym does not vanish in the 
past direction. This result is inconsistent with our initial assumption of f^YM = 0- Hence we conclude that the Kasner 
circle is not a real past attractor, in other words, there is no asymptotic solution that approaches the Kasner circle. 

We have confirmed this result numerically. In Fig. [SI we depict the time evolution of E± and the behavior of 
ab/H, bc/H, and ca/H. In this figure, we start with the initial data near the Kasner solution (point A: E = 0.81 and 
\f r = 0). Because Kasner spacetime has unstable modes, the spacetime does not approach the Kasner circle (E = 1), 
but rather E is going to decrease (Fig. |HJ (b)) and bc/H diverges (Fig. [S] (c)). Then the orbit again approaches 
another point of the Kasner circle. We may expect that the same process will be repeated because the Kasner circle 
is not an attractor. 

The mechanism preventing Oym — > is similar to that in the Mixmaster solution. In the case of the Mixmaster 
universe, when a spacetime approaching one point (ip) on the Kasner circle, one of N a grows while the others of N a 
damp. From the stability analysis of the Kasner circle in the vacuum Bianchi IX spacetime, the Kasner circle has 
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a one-dimensional unstable manifold, that is, N\ is unstable while N2,N^ are stable when — ^ir < ip < ^ir. When 
a spacetime approaches other points on the Kasner circle as ^ir < ip < tt, N2 is unstable, while N3 is unstable, if 
— 7r < tp < — . These are described by the following equations: 

Ni cx exp [2{2cos(V> + tt) + 1} r] 



N2 cx exp 



2 <^ 2 cos 



1 



3 cx exp 



2 ^ 2 cos 



yj + -7i ) +1 



(4.4) 
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FIG. 8: An orbit which starts near some point A on the Kasner circle evolves into another point B on the circle in the past 
direction: 

(a) The orbit in the E+-E_ plane. The initial values are E = 0.81 and \& = (point A), which is close to the Kasner circle (a 
solid circle). ~*S> eventually goes close to n (point B, which is slightly less than n). 

(b) The time evolution ^ for the same orbit. The orbit endures for a long time near the Kasner circle, and eventually it moves 
to another point on the Kasner circle. 

(c) The same orbit in the E+-(ab) 2 / H 2 , E+-(fec) 2 / H 2 , E+-(ca) 2 / H 2 planes. The initial values are E = 0.81 and (ab) 2 /H 2 = 
3.60 x 10 -7 , (bc) 2 /H 2 = 5.76 x 10~ 2 , (cafjH 2 = 1.60 x 10" 7 . * ~ results in the initial growth of (bc) 2 /H 2 , but * ~ tt 
(slightly smaller than tt) results in the decay of (be) 2 / H 2 and the growth of (ab) 2 /H 2 in the final stage. 



From Eqs. (|4.3(l and l]4.4[l . we see that the Bianchi I EYM model is quite similar to the Mixmaster universe if we 
replace Nx, N2, and with bc/H, ca/H, and ab/H. Therefore, although we cannot show the chaotic behavior in 
the EYM Bianchi I model in the past direction, we conclude that such a spacetime does not approach the Kasner 
solution but that it probably shows a chaotic behavior, just as the Mixmaster universe does. This indicates that 
energy density of the Yang-Mills field is not negligible compared to those of the shear and curvature terms. This 
breaks the second statement of BKL conjecture. 



B. Bianchi type IX spacetime 



Next, we consider the case of Bianchi IX spacetime. In this case we expect two types of chaos, i.e., the chaos of the 
Yang-Mills field and the chaos of spacetime. Do two types of chaos enhance each other, or evolve independently? A 
study of the Bianchi IX EYM model by Belinskii and Khalatnikov 12| gives hints on this question. They considered 
one Kasner epoch and found approximate solutions. Analyzing approximate solutions, it was shown that, although 
the spacetime is not affected by the Yang-Mills field between two Kasner epochs, the Yang-Mills field changes sharply 
in the transition. These results anticipate that the chaos of spacetime does not depend on the Yang-Mills field. In 
our study, it is very difficult to judge because numerical error causes our simulation to stop just after a spacetime 
approaches some point on the Kasner circle. Hence we shall take a guess by showing some indirect facts. Eq. I|4.1() 
yields Eq. (|4.3|l . when the spacetime approaches some point on the Kasner circle (Hym = 0). We also have 



(iVia) 2 cx exp[2{2 cos(V> + tt) + l}r] 



{N 2 b) 2 oc exp 



2 \ 2 cos j ip - ^tt 



I t 
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(N 3 c) 2 



oc exp 



2{ 2 cos ( -7T ) + 1 }t 



(4.5) 



and 



N 2 tx exp [4 {2 cos (ip + it) + 1} r] 



A 2 oc exp 



4 <j 2 COS ( V ~ gTT ) + 1 } T 



N% oc exp 



N1N2 oc exp 



1 <j 2 COS ( ( • • -- I 7 1 ) T 



I t 



4 ^ — COS ( tfj + — 7T 



iV 2 -/V 3 oc exp [4 {- cos(V> + 7r) + 1} 



N3N1 oc exp 



4 <j — COS ( tfj — -7T 



I t 



(4.6) 



These asymptotic behaviors show that SIym diverges as max [(A^a) 2 , (N 2 b) 2 , (N 3 c) 2 , (bc/H) 2 , (ca/H) 2 , (ab/H) 2 ] , 
while K diverges as max [Nf, N 2 , A/|, N1N2, N2N3, N3N1] . This means that even if a spacetime approaches one 
point on the Kasner circle once, it will leave the Kasner circle. We expect that it may approach another point on the 
Kasner circle, just as in the vacuum Bianchi IX spacetime, mainly because of the divergence of K. This is because the 
curvature term K diverges more strongly than the energy density of the Yang-Mills field f^YM • The growth rate of K 
is larger than that of Oym- The Yang-Mills field may not play any important role in this dynamics, which supports 
the results of Belinskii and Khalatnikov [12|. This means that the chaos of spacetime is stronger than that of the 
Yang-Mills field. This also means that the second statement of the BKL conjecture is satisfied, because SIym/K — > 
toward initial singularity. Considering that Bianchi IX spacetime is generic in spacelike homogeneous spacetime, we 
conclude that BKL conjecture is valid in the EYM system. 



V. CONCLUSIONS 



We studied the EYM system in class A Bianchi spacetime in order to g eneralize the result in axi-symmetric Bianchi 
I spacetime 0, ^] and that of the Einstein-Maxwell system 0, 0, |20| . We are also interested in the multiplicative 
effect of two types of chaos, that is, chaos found in Yang-Mills field and that found in Bianchi IX spacetime. We 
analyzed the case in which the shear and the Yang-Mills field can be diagonalized. We find that chaotic behavior 
appears in the late phase (the asymptotic future) for the Bianchi I, II, VLj, and VIIo- In this phase, the Yang-Mills 
field behaves as that in Minkowski spacetime, in which we can understand it by a potential picture. For types VIII and 
IX, a potential picture may no longer valid. Although we need further analyses, we find some sensitive dependence 
to initial data. It could be a new type of chaos in the Yang-Mills field. 

While, in the initial phase (near the initial singularity), we numerically find that it seems to approach the Kasner 
solution. However, we show that the Kasner circle is unstable and the Kasner solution is not an attractor. From the 
analysis of stability and our numerical simulation, we expect that it behaves as a Mixmaster universe even for Bianchi 
I spacetime. This result may provide a counterexample to the so-called BKL conjecture. However, it is not the case 
for the Bianchi IX spacetime, in which the BKL conjecture is still valid. 

We also analyze a multiplicative effect of two types of chaos. Two types of chaos seem to coexist in the initial 
phase. However, the effect due to the Yang-Mills field is much smaller than that by the curvature term. 

Finally we mention the transient phase. It is found that the behaviors of shear scalar X 2 and energy density of the 
Yang-Mills field Oym are very complicated. The potential picture of the Yang-Mills field is of course no longer valid. 
Hence this complexity is not only due to the chaos the Yang-Mills field appearing in the late phase. The chaos in the 
initial phase is understood by a mechanism similar to the Mixmaster universe. In the transient phase, however, the 
orbit does not need to approach the Kasner circle. So the behavior in the transient phase is not like the Mixmaster 
type. This could be new type of chaos that appears by a multiplicative effect, although it is a transient phenomena. 
This is under investigation. 
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APPENDIX A: DYNAMICS OF OFF-DIAGONAL COMPONENTS 

The vector potential of Yang-Mills field in a Bianchi spacetime is written as A = A a A ' (t)fA<^> a , using the gauge 
freedom. In this appendix, we show that if a ab , A a A ^ and Aa are initially diagonal, the off-diagonal components 
will not appear anytime in class A Bianchi spacetimes, while they will appear in class B Bianchi spacetimes. 

In class A Bianchi spacetime, each component of the Yang-Mills equation F .a — Ap x F a/3 = is written as 
follows: 

a = component : 

Ax x Ax + A 2 x A 2 + A 3 x A 3 - 2Q 3 Ax xA 2 - 2ft 1 A 2 x A 3 ~ 2ft 2 A 3 x A\ = 0; (Al) 

a = 1 component : 

Ax = -3HAx + 2Q 3 A 2 - 2fl 2 A 3 

+{-(<5- n + H) + (axx - 2H)(a u + H) + {a 12 + f> 3 ) 2 + (<j 13 - tt 2 ) 2 - (n 1 ) 2 }A 1 

+{{n 3 - <7i 2 ) - On - 2H)(a 12 - 3 ) + (a 12 + n 3 )(cr 22 + H) + (<r 13 - n 2 )(a 23 + n 1 )}A 2 

+{-(o 2 + 0-13) + On - 2F)(<7i3 + n 2 ) + O12 + n 3 )((j 23 - fi x ) + (a 13 - n 2 ){<r 33 + H)}A 3 

-(m + n 2 + n 3 )A 2 x A 3 - A 2 x (A x x A 2 ) - A 3 x (A x x 1 3 ); (A2) 

and a = 2[or 3] component is obtained by exchanging indices in Eq. I|A2I) as (1, 2, 3) — > (2, 3, l)[or (3, 1, 2)]. In class 
A Bianchi spacetime, some Jacobi identities for e a give a relation between Q a and a a t as 

n 2 + n 3 

III = 023 

n 2 - n 3 

n 3 + ni 

"2 = C31 

n 3 - ni 

il 3 = CT12 ■ (A3) 

ni - n 2 



The equation for the shear a ab is given by 

a ab = -3Ha ab + 2e cd {a a b)c n d - b ab + ^b c c 6 ab + ir ab} (A4) 

where 

b ab = 2n a c n cb - n c c n ab . (A5) 

We then assume a a b, and A a A ^ do not have initially off-diagonal components. From Eq. (|A3|) . Q a = 0. With 

this ansatz, Eq. (|A1() is trivial and Eq. I|A2I) is rewritten as : 

li = -3HA 1 + {-(& 11 +H) + (<T 11 -2H)(a 11 +H)-(n 1 ) 2 }A 1 

-(ni +n 2 + n 3 )4 2) 4 3) n - A^A^A^n - A^A^A^r,, (A6) 



Note that Q a — from Eq. I|A3(1 . The equations for A 2 and A 3 are obtained by exchanging indices. 
Setting a ab = for a ^ b and A a A%> = 0, 4 = for a =/= A, we obtain 

A^ = (ojtA). (A7) 

The off-diagonal components of dynamical equations for a ab is reduced to 

v a b = K ab (a ^ b) (A8) 
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because off-diagonal components of b ab vanish, where Tr a b is the traceless part of the energy momentum tensor of the 
Yang-Mills field, i.e., 7r ab = T ab — ^T^i] ab with T a/3 = F ai F^ — \g a pF 2 . Writing down the off-diagonal part of Ti ab 
explicitly, we can easily show it ab = as follows: 

7Ti2 = Fi 1 F 2 1 = F 10 F 2 a + Fi 3 -F 2 3 

= + (a n + H)A X ] ■ [A 2 + (a 22 + H)A 2 ] + [n 2 A 2 + A^A ( ^r 2 ] ■ [-n^ - A^A^f,} 

= 0. (A9) 

We also find ir 23 = and tt^i = 0. Thus we obtain 

& ab = (a^b). (A10) 

Eqs. (|A7|) and (|A10|) guarantee that the off-diagonal components of a ab , Aa and A^ always vanish. We conclude 
that there is some set of initial data in class A Bianchi spacetime for which a ab , A^ and A^ are always diagonal. 

Next we discuss the dynamics of off-diagonal components in class B Bianchi spacetimes. For simplicity, we consider 
only the Bianchi V spacetime, in which n a = and a a = (a, 0, 0). The basic equations for Yang-Mills field are given 
as follows: 
a = component 

2a(A 1 - n 3 A 2 + n 2 A 3 + (a n + H)A X + <r 12 A 2 + CT13A3) - Ai x A x - A 2 x A 2 - A 3 x A 3 

+2fi 3 i*i x A 2 + 2QxA 2 x A 3 + 2n 2 A 3 xl 1 = 0; (All) 



and a — 1 component is given by Eq. i|A2fl with n a = 0. 

a = 2 [or 3] component is derived from the corresponding equation in class A Bianchi spacetimes by inserting n a — 
and by adding the term 



Some Jacobi identities yield 



The equations for shear a a b is 



a 2 A 2[m . 3] + 2aAx x A 2[m 3] . (A12) 



a = -(H + a n )a (A13) 

n 2 = (73i (Ai4) 

n 3 = -a 12 . (A15) 



b ab = -3Ha ab + 2e cd {a a b)c n d + ir ab , (A16) 



where we assume a ab , Ai A) and Ai A) are diagonal. From Jacobi identities, we find Q 2 — Q 3 = 0, but we do not know 
about fii. We then rewrite Eq. (|A11|) as 

n x (A 2 x A 3 ) + a [A x + (<t u + h)A x ] = . (Ai7) 

Because the second term of the r.h.s. in this equation does not vanish (unless A^ = 0), we find fli 7^ 0, which 
results in ir 23 7^ 0. Hence the source term of the off-diagonal components of Eq. I|A16(I does not vanish and then a ab 

(31 ' (31 (2) ' (21 

becomes no longer diagonal. The Yang-Mills equations show that A\ , A 2 , A 3 , and A 3 ' become also non-trivial 

due to the term (|A12|) . We conclude that in type V Bianchi spacetime we cannot assume that a a b, A^ and A^ are 
diagonal even though they do not have off-diagonal components initially. 
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